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Introduction. — The theory of the incidence of waves of sound in a non- viscous air 
upon small obstacles of cylindrical or spherical form is well known to students of 
mathematical physics ; it has been treated in Lord Rayleigh's ' Theory of Sound/ 
and in Prof. Lamb's * Treatise on Hydrodynamics. 7 The corresponding problems for a 
viscous air have not, however, been worked out, and this paper is devoted to an 
investigation of these problems. The solutions of the equations of vibration of a 
viscous gas with reference to cylindrical and spherical surfaces were given by 
Prof. Lamb in a paper entitled " On the Motion of a Viscous Fluid Contained in a 
Spherical Vessel " and published in the ' Proceedings of the London Mathematical 
Society' in 1884. It is easy to obtain solutions suitable to the case of divergent 
waves ; the functions involved are Bessel functions with a complex argument. An 
analytical expression for the secondary waves diverging from the obstacle is obtained 
without difficulty. It then remains to find an expression for the loss of energy to the 
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primary waves. In calculating this loss of energy it is necessary to consider the 
dissipation of energy by friction in the immediate neighbourhood of the obstacle in 
addition to the energy which is carried away to a distance by the secondary waves. 
This was pointed out to me by Prof. Lamb, at whose suggestion this paper was 
written. In obtaining an expression for the energy dissipated by friction I at first 
made use of the dissipation function. This method led to exactly the same results as 
that finally adopted, but the mathematics involved were cumbrous, and the physical 
ideas, on which they were based, were somewhat obscure. Another disadvantage of 
this method was that it was necessary to calculate separately the scattered and the 
dissipated energy. I have to thank Prof. Lamb for his kindness in pointing out to 
me the method of calculating the lost energy adopted in this paper. The result has 
been to make the paper more clear and readable. 

I have succeeded in obtaining expressions for the energy lost to the primary waves 
in the case of spherical and cylindrical obstacles. As might be expected, the problem 
of the cylindrical obstacle presents greater analytical difficulty than that of the 
spherical obstacle, and in the former case it is necessary to obtain different 
approximate expressions according to the diameter of the obstacle. The results for 
wires of lO -1 cm. radius and for wires of 10~ 3 cm. radius can be obtained without 
much difficulty, but when the radius of the wire is of order 10~ 2 cm. it is necessary 
to perform very laborious calculations in order to arrive at intelligible results. The 
energy lost to the primary waves is, in all cases, very great compared with what 
would be lost in a non-viscous air, but the ratio of the lost energy to that incident 
upon the obstacle is at most of order 1CT 2 . 

In the case of spherical obstacles the difficulties of approximation are not so great, 
ss in the case of cylindrical obstacles the loss of energy is far greater than in a non- 
viscous air, but, as before, the ratio of the lost energy to that incident upon the 
obstacle is at most of order 10~ 2 . 

It is possible to extend the results obtained for a single obstacle to the case when 
the waves of sound are incident upon a large number of similar obstacles. This has 
been done by Lord Eayleigh for the corresponding problem in a non- viscous air ; the 
same method has been adopted in this paper. It should, however, be borne in mind 
that the results so obtained are valid only when the obstacles are so sparsely 
distributed that the space occupied by the obstacles is a small fraction of the total 
volume. The investigation has some practical interest. It has been asserted that 
the suspension of a large number of parallel wires in a hall or lecture room will 
improve the acoustic properties of the room. According to the ordinary theory of a 
non- viscous air the effect of any such arrangement of wires would be inappreciable. 
From the results of this paper it also appears that the viscosity of the air is not 
sufficient to account for the alleged phenomenon. 

The results in the case of spherical obstacles are of greater interest, since they are 
applicable to the consideration of the effect of foggy weather upon the propagation 
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and audibility of sound. If the diameter of the drops of water in a dense fog is 
assumed to be *02 mm., there is no appreciable alteration in the audibility of sound ; 
but, if the diameter of the drops of water is *002 mm., the presence of fog is distinctly 
prejudicial to the audibility of sound. The former case is in agreement with 
Tynd all's observations on the subject. 

In conclusion I desire to thank Prof. Lamb for much kind advice and encourage- 
ment in the writing of this paper. 

§ 1. In a viscous gas, if u> v 9 w be the components of the velocity at any point 
x, y> z of the fluid referred to fixed rectangular axes, and if p be the pressure at this 
point, the equations of vibration may be written in the form 



du 
~dt 

dv 



1 dp 

p dx 
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where p is the equilibrium density, v is a small constant of dimensions L 2 T -1 
depending on the viscosity, 3- has been written for div (u, i\ w)> and terms of the 
order of the square of the velocity have been neglected. 

By a method very similar to that # used in the case of an incompressible fluid it is 
found that the total rate of dissipation of energy within any closed surface S is 
given by 
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where q is the resultant velocity at any point of the fluid, Sn denotes an element of 
the normal to the surface S, l 9 m, n are the direction cosines of this normal drawn 
inwards in each case from the surface element 8S. Further, a, b, c, £ r), £ are given 
by the relations 
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When there is no motion of the gas parallel to the axis of z, and the motion is the 
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same in all planes perpendicular to this direction, the expression (2) for the dissipation 
of energy per unit length of the axis of z takes the form 



2F = 4 



3*V>0 



{a + b) 2 dx dy + ipp I £ 2 dx dy 



d(f 



+ 2vp (lu -f mv) (a + b)ds— vp \ y^- ds + &vp (Iv — mu) t>ds . . . (4) 

where ds is an element of the curve bounding the region in question. 

§ 2. We now proceed to obtain a solution of the equations of motion which shall be 
applicable to the case when the motion in all planes perpendicular to the axis of z is 
the same, and when further there is no motion parallel to this axis. 

In this case the equations of motion take the form 
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where 5- = — + — and V t 2 denotes the operator ^—^-h^—. 

dx dy ox ay' 

The equation of continuity takes the form 
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dt 
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(2), 



where squares of the velocity and other quantities of the same order have been 
ignored, and s denotes the condensation. 

If we neglect the effects of conduction and radiation of heat, we may write 



P = Po + p S 



( oK 



where p is the equilibrium pressure and c is the velocity of sound. 

Eliminating p and s from equations (1) with the help of (2) and (3), we obtain 
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These equations will be satisfied by 



u = 3 i + ^ 
9cc 3y 
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(5), 
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provided <f> and \jj are functions satisfying the equations 

l^cV^+lviv^ (6), 



a^ 

u 



vV^xjf (7). 



If we assume a time factor e t(rt , these equations take the forms 



(V^+Zi^^O, (V 1 2 +F)^ = (8), 

where h 2 and k 2 are given by 

h 2 = o-7(c 2 + |^o-), tf = -urfv . . . . . (9), (10). 

We shall for convenience suppose & to be equal to that root of equation (9) which 
reduces to cr/c when v is zero ; h will be taken to be equal to (cr/^) 1/2 .e"™ t7r/4 . 

With these conventions the solution of the equations (8), which represents waves 
of sound diverging from the origin, is given by 



00 ~\ 



<£ = A D (hr) + S A. n D n (hr) cos (n3 +' <% n ) 



00 



r* • • * t » I ill} 



xjj = B D (/cr) + S B B D B (&r) sin (n£ + /J n ) 



n~l 



where for convenience the time factor has been omitted and where D n (£) is given by 

J n (^) and Y n (£) have their usual significance as Bessel functions, and A n , B m a n , and 
/3 n are arbitrary constants. 

It need hardly be remarked that to obtain the actual expressions for <f> and \fi it is 
necessary to multiply the expressions contained in (11) by e la \ and to equate <j> and \}t 
to the real parts only of these products. For the sake of brevity we shall usually 
omit, when possible, the time factor e l<r \ 

Since, in the case of air, v is a small quantity, when expressed in cm. sec. units, it 
is clear from (9) and (10) that for all audible sounds \h\ will be large compared with 
h\. In fact, at any ordinary distance r from the origin the \fi terms in (11) will 
become insensible owing to the factor exp. { — (^a m /v) 1/2 r}. 

For plane waves propagated in the negative direction of the axis of x 9 the solution 
will be given by 

<£ = 0*V, xjs = ........ (12), 

where, as before, only the real part of <f> is to be taken into account. 
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§ 3. The Incidence of Plane Waves of Sound upon an Obstructing Cylinder. — We 
are now in a position to consider the effect of a cylindrical obstacle upon a, train of 
waves propagated in a direction perpendicular to the axis of the obstacle and incident 
upon it, 

We take the axis of the obstructing cylinder as axis of 2, and suppose the incident 
sound to be propagated in the negative direction along the axis of a?. Then, as in (12) 
of the last article, we may assume for the incident sound the expressions 

<£o = e** t </r = 0. 
Expanding in series of Bessel functions, we obtain 



00 



(j) = J (hr) + 2 2t n J n (hr) cos nS-, i// = . . . . . (1). 



n = l 



The scattered waves will be symmetrical about the axis of x 9 or £ = 0, and 
consequently we may assume for them the forms 



co 



fa = A D (hr) + % {k n "D n (hr) cos nS }, 



n- = 1 

QQ 



fa = 2 {B M D tt {hr) sin n&). (2). 



n = 1 



At the surface of the obstacle the radial and tangential components of the velocity 
must vanish ; hence we must have 

when r = a, if a is the radius of the cylinder. 

In order that the boundary conditions (3) may be satisfied, we must have 

AJiaDJ (ha) = —ha J J (ha) 
or 

A D 1 (^a) = — J 1 (ha) (4), 

and in general for n>0 



AJia DJ (ha)+nB H D n (ka) = —2i n ha3J (ha) 
nk. n D n (ha) + BJca D/ (to) = — 2u n nJ n (ha) 



• (5). 



These equations (4) and (5) are sufficient to determine the various constants in the 
expressions (2) for the scattered sound. In the process of approximating to the values 
of these constants by means of equations (4) and (5), we shall confine ourselves to the 
case when ha is a small quantity ; in other words, we shall assume that the dimensions 
of the obstacle are small compared with the wave-length of the incident sound. The 
other case, when the dimensions of the obstacle are large in comparison with the 
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wave-length, presents exactly the same difficulties as occur in the similar problem in 
connection with the theory of a non-viscous gas. 

We shall also consider especially the case when the viscous gas is the air of the 
atmosphere; in this case v is a small quantity about '132 in cm. sec. units, and 
consequently we may regard cvjc 2 as a small quantity for all wave-lengths. Since 
av\(? is small, we may write very approximately from (9), § 2, 



■h = ^(1-funVc 2 ) 



(6). 



We must now obtain approximations to the values of the constants in the expressions 
for the scattered sound by means of equations (4) and (5). For this we shall need 
the approximate values of the Bessel functions involved in these equations. For 
convenience we shall write them down. 

When £ is small, we have for all values of ^>0, 

D n (£) = — '-^ £~ w + less important terms 

B n '(i)= -2^L ! t ( " +1) +... 

IT 

J n (£) == — — : £ n -f terms containing higher powers of £ 

J J (£) = --— _ ~ £ w-1 + terms containing higher powers of £ 

Li . (71/ 1 ) I 
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■\ 
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= —2i n nha { J n (ha) DJ (ha)—D n (ha) J n ' (ha)}. 



Further 



j. (o r>: (Q-D n (o j.' (4) = - 1 y„ a) j,' (o-yj (o j n (oi = - - r 

7T ' 7T 

by a well known result in the theory of Bessel functions. 
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Closer approximations are 

d.(£) = |{-(iogK+r+»+K 2 (iogK+r+^-l)} 

Also we have 

Jo(0 = i-K 2 , J 3 (0 = K(i-K 2 ), J 2 (0 = K 2 (i-An • • (9). 

Now, on eliminating A n from equations (5), we obtain 
B n \Jia . ka D n ; (/*a) D n ; (ka)—n 2 D n (ha) D n (ka)~\ 
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Consequently we obtain 

B» [ha . ha D,/ (ha) D,/ (ka) - n 2 D n (ha) D n (ha)] = — l" . . . (10). 

77 

Now ha is a small quantity, and consequently we may use the approximations (7) ; 
substituting for D n (ha) and D/ (ha) from (7), we obtain 

-^h^a^BJJcaBJ (ka) + nB n (ka)} = ~t w 

7T \ . ^ 

or 

_: t h~ n a~ n . ha D^ (&a) B n = — t n . 

7T 7T 

Hence we obtain as a first approximation 

■R -^ 1 h*a* ( . 

n ~~ 2 n ~ 1 ' (n-l)l' kaD n ^(ka) [ >' 

Similarly, by elimination of B n between equations (5), we obtain 

A n [ha . ha D' n (ha) D' n (ka) —n 2 D n (ha) D n (ha)] 

= — 2i n [ha . kaJ f n (ha) D' n (ka) —n 2 J n (ha) D n (ha)]. 

Using the approximate values given in (7), we obtain 

2 n n\ _ _ h n a n 

-A n — ^ n a n kaD n ^(ha) = 2t n ~— — — — D n+1 (&a) 
7T 2 n (n — 1)! 

or 

A _ _9," 2Z_! ^ J_Jj±l v 1 ^/ /19\ 

n ~ •2 2 ^!(r,-l)!"D B _ 1 (/^) ^ ' 

for all integral values of n > 0. 

By a similar process, but carrying the calculation to a higher degree of approxi- 
mation by means of (8) and (9), we have 

a = _i._ A 2^2 D 2 (ha) -jh 2 a 2 {2J) 2 (ha) -D (ka)} , , 

2 D (&a)-|AW(log|^ - { >' 

Let us consider first the case when | ka | is small. In this case it will be sufficient 
to derive the value of Ax from (12). Hence we have approximately 

z D (ka) 

Writing for convenience 

k = Xe~ 1/4l7r , where X = (or/^) 1/2 (14), 
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and using the approximations for D 2 (£) and D (£) given in (8), we have 



D u (ka) 



7T 



2 



(log \\a + y + £«r) - |-iW (log %\a + y + \ «r - £)] 



D 2 (*a) =-[2i\- 3 a- 2 +J]. 

Substituting these expressions for D (&a) and D 2 (ka) in the above formula for A 1; 
we obtain 



A i 



IT 



h 2 d 



1-iiW 



X 2 a 2 * log |Xa + y + itvr + i^X V (log |Xa -f 7 + %nr - 1-) * 



In general, it will tend to sufficient accuracy if we take 



A, 



(TV 



- T 7r{log^\a + 7 + ii7r} * 



• « 



c 



. . (15) 



in the case when Xa is small. 

Similarly we obtain from (11) the approximation 

Bl = ^^(^gi^+r+i^)" 1 

when Xa is small. 

Let us next consider the case when Xa is large. 
Since J Jca | is great, we may write 



» « * » 



. (16), 
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Substituting these expressions in the formula given for A x in (13), we obtain 
approximately 

f I61 112i 2 1 

A x = iiirh'a? { 1 - gfo + p^js -i^a 2 (log £fta + y+b^+i) j * 



When Xae 1;4t7r is written for jfca, this reduces to 

A x = ^nrh 2 a 2 [1 + ^/2 . (Xa)" 1 -\h 2 a 2 (log ^a+7-f f) 

-i{ v /2.(Xa)- 1 +l(Xa)" a + iirA 2 a 8 }] . . (17). 



As a first approximation to the value of B l9 we find in the case when Xa is large 



B x = 2i/ia^--|-Y / % t( ^ +1 ^ ) . 



« • • • 



(18). 



The approximate value for A , obtained from equation (4) with the help of (8) and 
(9), is easily seen to be 

A = -^7r/^a 2 {l-f|-A 2 a 2 (log|-/ia + 7 + 4-i7r— f)} . .... (19). 
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§ 4. Now that we have approximated to the values of the various constants involved 
in the expressions for the secondary waves, we can proceed to estimate the additional 
rate at which energy is being dissipated in the space surrounding the obstacle. This 
additional rate of dissipation will be equivalent to the rate at which energy is being 
lost to the primary waves in consequence of the presence of the obstacle. Now, if we 
consider a region bounded internally by the obstacle and externally by a cylindrical 
surface coaxial with the obstacle and of radius R, it is clear that the difference of the 
rates at which energy is being carried across the internal and external boundaries of 
this region will be equivalent to the rate of dissipation of energy within it. 

Hence, if p^ q denote the pressure and the radial velocity at any point due to the 
primary waves alone, and if p u q 1 denote the pressure and radial velocity due to the 
secondary waves alone, it is easily seen that the dissipation of energy within a 
distance R of the obstacle is given by 

(Po+Pi)(qo + qi)ds . . . (1), 

where it has been assumed that q , q 1 are both measured inwards, and the integration 
is to be taken round the boundary of the surface r = R. 

Now the dissipation of energy due to the primary waves alone is given by J jVZo ds. 

Hence it follows from (1) that the additional dissipation of energy due to the 
presence of the obstacle is given by 

p qids+ piq ds.+ 2 3 ^ds (2). 

Now the rate at which energy is being carried across the surface r = R is —jpiqi ds, 
and hence from (2) it follows that the total additional dissipation of energy due to 
the presence of the obstacle is expressed by 

p q 1 ds + \piq ds (3), 



4k 



where the integration is to be taken round the boundary of the surface r = R. # 

Since crvfc 2 is in all cases a small fraction, it is clear from (6), § 3, that R may be 
great compared with the wave-length of the incident sound and yet such that orVR/c 3 
is a small fraction. In this case we may neglect the imaginary part of JiR in 
expressing the value of <f> and <£ x at the surface r = R. Also, if crR/c is great, it 
follows that | Ml | or XR is very great, since their ratio is the small quantity (crv/c 2 ) 112 . 
Since, when | Ml | is great, J) n (kB,) approximates to the value 



JLY /2 

-Ml/ 



.n -i (AR/ V 2 + V 4 t) o ~ AI V ^ 2 
,77/ 5> ' 



* This method of finding an expression for the loss of energy was kindly suggested to me by 
Prof. Lamb. I had previously obtained the same results by means of the dissipation function ; but the 
work involved was very cumbrous. 
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we see that the \js terms in the expression for the secondary waves are inappreciable 
at the external boundary owing to the exponential factor exp. (—\R/ v /2). 

"We shall suppose then that R is great compared with the wave-length of the 
incident sound and yet such that <xVR/c 3 is small. 

Let us now return to the consideration of (3). At the external boundary, r = R, 
we may write approximately 

/ O \l/2 oo 

T \7rJlRJ „ = o L J 

where square brackets are used to denote that only the real part of the expression so 
enclosed is to be considered. 
Hence we obtain 

Pl = - po Ep= -po<r( -fp t [A n t n+1 cos n§ . e^-AB-^ir)] ? 

Ot \7rllXxJ n = 

^J. / 9 Xl/2 00 

qi = - Sp = A /jL $ [A^* 1 cos n£ . ^" M ™ 
Using these expressions for p x and ^ we obtain 



1/2 00 



M+.Mi = -(-Td) (po<Tqo-hp ) $ [A n L n+l cos n^e'^-* 11 - 1 ^)]. . . (4). 

Again, since AR is great, we have from (1), § 3, approximately 

/ 2 \ 1/2 . °° / 2 \ 1/2 

<£o = — To ) sin ( AR + -J^r) cos cr£ + 2 £ -y^ sin (KR +far— %nir) cos (cr£ + inir) cos wS. 

Hence we obtain 

/ 2 \ 1/2 . . °° / 2 \ 1/2 

p = /o cr —jHg- sin (AR + ^7r) sin cr£ + 2p cr 2 ( -7-5- ) sin (AR + ^Tr—^nrr) 

sin (<rf + |wr) cos nS- 

cLJJLvJ. 

/ 2 \ 1/2 °° / 2 \ 1/2 

g = —h( —j^ cos (AR + J?r) cos <rt--2h t -7-5- cos (AR+ J-7T— |mr) 

\7rh1iJ " »=i \7rhhjJ 

cos (cr£ + |rmr) cos n#. 
Combining these expressions for po and q we find without difficulty 

— (pQ<rqo—hp ) = hp a ( ~yp ) cos ((rt — KR—l7r) + 2hp (T (-r^ 



00 



S ( — ) w cos (cr^— AR— ^tt) cos nS-. 



n=l 



Substituting in (4) and integrating with respect to 5- we obtain 

*2rr 00 

(jPi2o+ j pog f i)RdA= 4p (rS [(~) n A n t tt+1 cos(cr^-AR-i7r)e l(<r< - AR - 1/ ^] 
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of which the mean value is 

00 

ZOqCF 2* (' — ) -A-ftk 1 ... . , , . V^/« 



n=0 



This last expression, then, represents the loss of energy to the primary waves in 
consequence of the presence of the obstacle. 

From the value of A n obtained in (12) of § 3 we see that the summation (5) consists 

of a series of terms arranged in descending order of magnitude. Consequently, in 

determining its value we shall limit our attention to the first two terms of the 

summation. Hence the total loss of energy to the primary waves is given very 

approximately by 

2p cr [t A + AJ (6). 

Let us first consider the case when \a is small. In this case we have, from (17) 
and (19) of § 3, 

A i = -^^(logi-Xa + y + liTr)"- 1 , 

A = - %nh 2 a 2 { 1 + \h 2 a 2 (log \lia + y + \nr - f ) } . 
Hence 

[ t A + AJ = 1 l7r%V-^7r(logi\a+ r ){(logi\a + y) 2 + 1 V 2 }- 1 - 

Now for small values of the radius a the first term of this last expression is small 
compared with the secondhand consequently may be neglected. Hence the loss of 
energy to the primary waves is given approximately by 

-2p ^- 2 7 r(logl\a+ r ){(logjX« + r ) 2 + 1 V 3 }- 1 (7), 

in the case when \a is a small fraction. 

The ratio of this last expression to p cr 2 a/c, which represents the rate at which 
energy is incident in the primary waves upon the obstacle, is given by 



^^(loglXa + r ){(loglXa + r ) 2 + 1 V 2 }" 1 (8). 



We may notice from this last result that, when \a is small, the proportion of the 
incident energy, which is lost to the primary waves, is very nearly proportional to 
the reciprocal of the radius, since the logarithmic terms will change more slowly as 
the radius changes. Hence, as in (7), the total energy lost to the primary waves is 
almost independent of the dimensions of the obstacle, provided these are small enough 
to satisfy the conditions under which the results (7) and (8) have been obtained, and 
provided also that we limit our attention to obstacles whose dimensions are of about 
the same order of magnitude. It might have been anticipated that the energy lost 
by friction in the neighbourhood of the obstacle would, in the case of very small 
obstacles, alter very slowly with the dimensions of the obstacle, and consequently 
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that the effect upon the primary waves of very small obstacles would be almost 
independent of the dimensions of the obstacles, provided they were limited to be of 
the same order of magnitude. 

Let us now turn to the case when Xa is great. In this case we have approximately, 
from (17) and (19) of § 3, 

A x = IittW [1 + ^/2 (Xa)" 1 -\h 2 a 2 (log %ha + y + \) ~l{^/2 (Xa)- 1 + f (Xa)~ 2 + J^a 2 }], 
A = - \Tih 2 a 2 { 1 + \h 2 a 2 (log \1\ a + y + \vn - f ) } . 

Hence we find approximately 

[1A0 + AJ = -^/iW^tt^^ 

- f^h'a" + l^a 2 { ^/2 (\a)~ l + 1 (Xa)" 2 } . 

Substituting this value for [iA -f- AJ in (6), we obtain, when Xa is great, for the 
total loss of energy to the primary waves the formula 



^pocrhW+Trpoatfa 2 {^/2 (Xa) ^f (Xa) 2 }. 



The ratio of this to p o- 2 afc, which represents the rate at which energy is incident 
upon the obstacle, is given by 

f^^aV^ + TTv^^V^c + l^ca) (9), 

which gives the proportion of the incident energy which is lost to the primary waves. 
The first term in (9) is independent of the viscosity of the medium, and is obtained 
in the ordinary theory of a non- viscous air. The second and third terms of (9) 
represent the additional loss of energy to the primary waves consequent upon the 
viscosity of the medium. Further, since the ratio of the second to the third term of 
(9) is of order Xa, it follows that the latter may be disregarded. Hence we see that, 
since a does not enter into the second term of (9), the additional proportional loss of 
energy consequent upon the viscosity is almost independent of the magnitude of the 
obstacles. In other words, the actual loss of energy in the primary waves due to 
friction is proportional to the radius of the obstructing cylinder if this be sufficiently 
large. This last result is clearly what might have been expected on physical 

grounds. 

It remains to consider the case when Xa is neither very small nor very great. In 
this case it is impossible to approximate to the values of the Bessel functions involved. 

From (13) § 3 we have 



A x = —\nrWa 2 



D 2 (ha) 



D (ka) ' 
Hence it follows that 

E 2 (Xae 1 ^) _ T> 2 (\ae- 1 '** ) 



[AJ = ^irh 2 a 2 . |t 



E (Xae 1 ^*) D (Xae" 1 ^) 

2 K 2 



> ? 
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where 



2 



E.(0 = -{fag 2- y +*«r) J„(0-Y„(£)}, 

7T 



and square brackets are, as before, used to denote that the real part only of the 

expression so enclosed is being considered. 

Hence we have 

[AJ = -i«rAV(D a Eo-D E a ) |D ( 



-2 



Now [iA ] = i- 6 -7r 2 A 4 a 4 , and consequently it may be neglected in comparison 
with [AJ. 

Hence the loss of energy to the primary waves is given by 

Jt7r/)o<rAV(DoE a -DJS )|D ( 



The ratio of this to p cr 2 a/c is 



h n 



aa 



(D E 2 -D 2 E )|D |- 2 (10), 



which therefore represents the proportion of the incident energy, which is lost to the 
primary waves. 

On p. 253 will be found a table giving the ratio of the lost energy to that incident 
upon the cylindrical obstacle in a number of different cases. For wires and 
cylindrical rods of comparatively large radius it is necessary to use the formula (9) ; 
the results for wires of radii 10 cm., 1 cm., and *1 cm. have been deduced from this 
formula. The formula (8) is applicable when the radius of the obstacle is very small, 
and the results for wires of radius 10~ 3 cm. have been obtained from it. When the 
radius of the wire is of order 10~ 2 cm., neither of these approximate formulae is 
applicable, and it becomes necessary to calculate the results directly from (10) ; the 
value of the ratio of the lost energy to that incident upon the wire has been worked 
out in this case for only a few values of the wave-length on account of the laborious 
nature of the work involved. 

It should be added that in the table given on p. 253 X denotes the wave-length of 
the incident sound measured in centimetres, and K denotes the ratio of the lost 
energy to that incident upon the obstacle. 

I have also worked out the values of K for wires of different diameters in the case 
when the wave-length of the incident sound is 250 cm, The results are arranged 
below : — 







A = 250 cm. 










1 cm. 

• 
• 


•1 cm. 

• 


•05 cm. 

• 
• 


•01 cm. 

• 
• 


•005 cm. 

• 
• 


•001 cm. 

• 
• 


K . . . . . 


• 
• 

•16 lO- 2 


• 
• 

•15 10- 2 


• 
* 

•21 10-2 


• 
• 

•26 10~ 2 


• 
• 

•32 lO- 2 


• 
• 

-8Q lO" 2 
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§ 5. Extension to the Problem of a Number of Cylindrical Obstacles. — Let us 
suppose that there are n parallel wires per unit area of a plane perpendicular to their 
common direction. At a distance r from the axis of any one of these wires, great 
compared with the wave-length of the incident sound, the secondary waves due to 
that wire will be given very approximately in all cases by 

<f) x = A D (hr) + A x D x (fir) cos $. 
Since hr is great, we may write for all values of n 



D, (hr) 



2 



\l/2 



IT 



hr, 



t n^-i(Ar + V*w) 



and ^>! takes the form 



(/>! - (A + A x t cos a) (^) 1/2 e -^ hr+l ^\ 



Along the course of the primary waves (3- = rr) this reduces to 



(!)• 



Consider now the wires which occupy a thin stratum dx perpendicular to the course 

of the primary waves. Let AP be the section of this 
stratum by a plane at right angles to the wires and let be 
the point, at which the vibration is to be estimated at a 
P great distance from the stratum. 

If AP = z, the element of area is dx . dz y and consequently 



the number of wires cutting it is n . dx . dz. 

Also, if OP = r 9 AO = — x, then r 
zdz = r dr. 

The resultant at of all the secondary vibrations, which 
issue from the stratum, is by (1) 




x 2 + z 2 and 



* cc 

2ndx j (Aq—Axi) 

+• —x 



/ O \l/2 



,7T 



h, 



g-i(Ar + V 4 »r) 



r V2 fa 



^/(r 2 —x 2 )' 



Writing r = —x + rj, we obtain 

2 



In evaluating this last expression we may assume rj/x to be very small ; it then 

takes the form 

2 X1/2 



Ziiv ax 



IT 



h. 



(Ao- A lt ) e'<**-v.-> (2 v )- 1J2 e-^ d v . . 

Jo 



» « 



• (2). 



Now 



-CO /»CO 

Jo Jo 
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Also by a well known result we have 

r°° r 00 / —. \i/2 

cos (hv 2 ) dv = sin (hv 2 ) dv = ( — ) 

Jo v .0 ' \ollJ 

Consequently 



JW \ -1- / W 



Substituting this result in (2) and restoring the time factor, we obtain for the 
resultant at of all the secondary vibrations coming from the stratum dx 

-2n.dx.(A 1 + LA )h~ 1 e l{hx+crt) (3). 

When Xa is small, A x is great compared with A ; neglecting A and using the 
expression for Ax given in (15), § 3, we obtain instead of (3) 



vrr 



2n .dx.-f (log £Xa+y+£Mr)~y (Aar+,rf) , 



of which the real part is 



PIT 



2n dx — { (log -|Xa + y) cos (hx + crt) + \tt sin (hx + at)'} / { (log jXa + yf + yqTt 2 } (4). 

To this is to be added the corresponding expression for the primary wave 

<£ = cos (Jix + a-t). 

The coefficient of cos (hx + at) is thus altered by the obstacles in the layer dx from 
unity to 

1 + 2n dx ~ (log !Xa + y)/{ (log %ka + y) 2 + -^tt 2 } . 
Thus, if E be the energy in the incident waves, we have 

cffi/E = ±n.dx — (log |Xa + r ) /{(log %\a+yf+-fo?}. 

Integrating this, we obtain 

E = E e~~ ax 9 

where E is the energy in the primary waves at incidence and a is given by 

a ^-i^f(\og^a + 7 )/ {(log iXa + yf+^W} (5)- 

The coefficient of sin (hx + <rt) in (4) gives the refractivity of the medium as 
modified by the wires. If S be the retardation due to the wires of the stratum dx, 

3 = -h^dxl {(log iXa+y^-HV 2 }. 
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Hence, if /x be the refractive index as modified by the wires 



Hence we have 

/x 



^-i = i^ 2 ^/{(iogi^+r) 2 +iV 2 } (6). 

1 = i- 7r ^^/{( 1 °gi Xa+ r) 2 +iW} (7) 

era J 



where p denotes the ratio, assumed small, of the volume occupied by the wires to 
the total volume. 

Let us now consider the case when Xa is great. In this case we have, from (17) 
and (18) § 3 

Aj + iAo = \nrh?a 2 [1+2^/2 (Xa)~ 1 --|A 2 a 2 (log %h,a + y+-£%) 

-i {2^/2 (\a)- 1 + i(\a)~ 2 +frh 2 a 2 }']. 

Substituting this expression for Ax + iAq in (3), we obtain for the resultant at of 
all the secondary vibrations coming from the stratum dx 

fa = —^inir dx . ha 2 [1 + 2 ^/2 (Xa)" 1 — f h 2 a 2 (log \ha + y+ ^) 

-i {2 <J2 (XaT + f (Xa)- 2 +f ttAV}] &<**+** 
of which the real part is 

— \mr . dx . ha 2 {2^/2 (Xa)~ J + f (Xa)~ 2 + f7r^ 2 a 2 } cos (hx + crt) 

+ ^mr.dx.ha 2 {1+2 ^/2 (Xa)~ 1 —^h 2 a 2 (log^ha + y + ~-^)} sin (hx + crt) . (8). 

To this is to be added the expression for the primary waves 

<j) = cos (hx + crt). 

The coefficient of cos (hx + crt) is thus altered by the obstacles in the layer dx from 
unity to 

l—\mrha 2 {2 </2 (\a)~ l + % (Xa)~ 2 + §7rA 2 a 2 } dx. 

Hence, if E be the energy in the incident waves, we have 

<£E/E = -mrha 2 {2^2 (\a)~ l + l (\a)" 2 + ^Trh 2 a 2 } dx. 

Integrating this, we obtain 

E = E e-*, 
where 

a = nnha 2 {2 ^2 (Xa)" x + f (Xa)~ 2 + f wh 2 a 2 }. 

When cr/c and (cr/*>) 1/2 are substituted for h and X respectively, this takes the form 

a = 2 ^72^77 . a . cr 1/2 v 1/2 /c + Imrvjc + f n7r 2 cr 3 a 4 /c 4 ( 9 ). 
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The second term of (8) gives the refractivity of the medium as modified by the 
wires. If S be the retardation due to the wires of the stratum dx 



8 = \wna 2 < 



1 + 2 ^2 (Xar-f^log *o»/c+ y + A) 



Hence, if /x be the refractive index as modified by the wires, we have approximately 

p-l = ^{1 + 2^2 (\a)- l }p (10) 

where p denotes the ratio, assumed small as before, of the volume occupied by the 
wires to the total volume. 

If the waves of sound traverse a medium in which a number of parallel wires are 
arranged, then the reciprocal of a will determine the distance which the waves will 
travel before the intensity of the sound is diminished in the ratio of 1/e. For sound 
waves of wave-length 10 cm. passing through a medium, in which there are 100 
parallel wires of radius 10~ 2 cm. per unit area of a section perpendicular to the wires, 
this distance is 47 cm. For greater wave-lengths the distance is greater. It seems 
hardly probable that any arrangement of wires could improve the acoustic properties 
of a room unless some other factor than viscosity is taken into account. Of course, if 
n is made sufficiently great, the reciprocal of a may become very small ; but it seems 
probable that it would be difficult to arrange the wires so closely that n should be 
greater than 10 3 . If it was possible to arrange wires of radius 10~ 3 cm. so closely that n 
was 10 4 , then the intensity of sound of wave-length 40 cm. would be diminished in 
the ratio 1/e after passing through a thickness of less than 4 cm. Such a contrivance 
could hardly, however, be carried out in practice. 

§ 6. Problems Relating to Spherical Obstacles. — We require first a solution of the 
equations of motion suitable to such problems. Differentiating the equations of 
vibration (1), § 1, with regard to x, y, z respectively and adding, we obtain with the 
help of (2), § 2, y a 

ar cVs+3 " v dt { h 

If we now assume a time factor e l<r \ this equation takes the form 

(V 2 + h 2 )s = 0, where h 2 = o- 2 /(c 2 + f^cr) .... (2), (3). 

Also the equations of motion (1). § 1, may with the help of (2) and (3), § 2, be 
written in the form 

(V*+1^u = (W-h*)Q, ^+k 2 )v = (k 2 -h 2 ) d ^-, (V 2 +k 2 )w = (k 2 --h 2 )^ . ( 4 ), 

where 

(j) = icrh~ 2 s = —h~ 2 div (u, v 9 w) and k 2 = —icr/v . , . (5), (6). 

These equations (4) are satisfied by 

(u t v, w) = grad <f>, 

VOL. COX. — A. 2 L 
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where <f) is any solution of the equation 



(V 2 +A 2 )(£ = (7). 



The complete solution of the equations of motion will be given by 

(u y v, w) — (u\ v\ iv f ) + grad <j> . . , 



(8), 



where u\ v\ w f satisfy the equations 



together with 



(V 2 + k 2 ) u f = 0, (V 2 + If) v f = 0, ( V 2 + h 2 ) iv f = 



du f dv f dlV f __ n 
dx dy dz 



/* « » » \ *j I* 



The solution of these equations suitable to the case of waves diverging from the 
origin is given by 



u' = $\(n+ 1)/,-! (kr)^~nk 2 r 2n +% +1 (Jcr) 



a 



O) 



n 



fi x \ r 2n + l 



, + 2/.(*r)(y|-*|.) x . . (10), 



with corresponding expressions for v f and iv\ Here o> M and x?i are solid harmonics of 
positive degree n, and j£ (hr) is a function of hr given by the relation 



MQ = [-~^ e ^- = %AQ-^{i) ■ 



• • 



• (ii). 



where *&„ (£) and \\i n (£) satisfy the relations 



*.(£) = ~ 



cl \ n cos £ 



4 ^b/ b 



^n (0 = ( ~ 



d \ n sin £ 



w . 



• • » » 



(12). 



The general formulae for «//„(£), ^ n (£). and /„(£), are 



$» (0 = 1 o e 



c 



2 



+ 



r 



(2w+l)l 2(2^ + 3) 2.4.(2» + 3)(2n+5 



:)-•)■ ' ' (13) ' 



T / y\ X . O • O ■ • • ( J-i ih ~~ " J- J 



1- 



£ 






? 



2(l-2n) 2.4.(l-2w)(3-2n) 



• . 



■ (14), 



/- (0 = 



i"e~ lf 



I 



» + l 



< 



1 + 



r^ 



. (?i+l) (n— 1) n(n+l) (n + 2) 



2i£ 






2.4.(i£) 2 



-f ... 



1. . Zi . O . . . LiTb 

2.4.6...2n(i£) 



ft 



• (15). 



The functions i|>„ (£), ^n (£)> an( i/» (£)> a ^ satisfy recurrence formulae of the types 

^U)=-&K + i(£)> 6/^(0 + (2w+ 1)^.(0 = *.-i(9 • (16). (17); 



these will be found useful hereafter in reductions. 
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Returning now to the consideration of the equations of vibration, we find from (8) 
and (10) that the general solution suitable to divergent waves is given by 






(n + l)f n ^(kr)^-nFr^% +1 (kr) 3 ( ^ 



00 






+ ?J^H-%)^ ^ 



with corresponding expressions for v and w ; <£ here represents the general solution of 
(7) suitable to divergent waves. 
Hence 

00 

<£ = t f n (hr).<f> n (19), 



tt=0 



where <\> n is a solid harmonic of positive degree n. 

If the motion is in planes through the axis of x 9 and is symmetrical about that 
axis, the solution takes the form 



u — -3- + 2 

OX 71=0 



'(n + l)f^(kr)^-nF^f n+1 (kr)^^)} . . (20), 



with corresponding expressions for v and w. Further we have 



00 



<£= t A n / M (Ar) r n P B ( M ), a, = B n r"P„ (^) . . . (21), (22), 



n = 



where A n and B n are arbitrary constants. 

We write as before 

h = Xe- 1 /^, where X = (<t» 1/2 .... (23), (24). 



Va 



From (15) it is seen thatj^&r) contains the exponential factor e~ {1Mv) %T or e~ 1/2A/2kr ; 
consequently since Xr becomes very great within a short distance of the origin, it is 
clear that, at a moderate distance from the origin, those parts of the expression (20), 
which depend on the functions f n (/cr), become inappreciable and may be neglected. 
Hence, at a sufficient distance from the origin, we may write 



co 



(u, v 9 w) = grad <£, where <j> = % A n f n (hr) r n V n (jjl). 



n=Q 



§ 7. The Incidence of Plane Waves of Sound upon an Obstructing Sphere.— We 
may now consider the effect of a spherical obstacle upon a train of plane waves of 
sound. Suppose the centre of the obstructing sphere to be at the origin, and the 
sound to be propagated in the negative direction along the axis of x ; then, as before, 
we may assume for the incident waves 

(u 09 v , iv ) = grad <£ , where <£ = e Lhx . . . . . (1), (2). 

2 L 2 
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Expanding in terms of the functions \fj n , we obtain 



00 



<£ = 2 (2n+ 1) i K /A//„ (hr) r n P n (fi) 

n=Q 



(3). 



The scattered waves will be symmetrical about the axis of x f and so we may assume 
for them 



w-0 



00 

(u u v u Wj) = grad ^+ S i (n+ l)f n -i (for) grad o) n —nh 2 r 2n ~ d f n+1 (hr) grad 



0) 



n 



X 



2n + l 



• (4 



where 



00 

<£, = 2 KJn (hr) r»P„ (/*), a> n = B B r»P B (/*) .... (5), (6). 



M = 



Now, by means of the recurrence formulae, it may be proved without difficulty that 



S- { f. (hr) r"~P n \ = — !— 



/.-i (M^P.) +AV> + »/. +1 (^)^(^i 



• (7). 



At the surface of the spherical obstacle r = a, we must have 



^0 + % = 0, Vq + ^x = 0, 



Wq + w x = 0. 



Hence, when r = a, we have 



n=0 



00 

grad(<£ + <£i) + 2 \(n+l)f n ^ (hr) grad co n ~-nh 2 7> 2n+ % +i (hr) grad 



0) 



tt 



2w + l 



,f 



= 0. 



Introducing the expressions given above for <£ , <jS x , to,,, we obtain for all values of n, 
when r = a, 

grad {(2ra+ 1) L n h n xh n (hr) r n V n ) + grad {A*/ B (/ir) r n V n } 

+ (n + !)/,_! (*r) B n grad (r*P n ) -ntfr* + *f n+l (hr) B n grad f ^ = 0. 

Hence, by means of the identity (7), we find that the following expression 



?h» k_> (ha) A (^p.) + /,v-% +1 M ^ (3 

+ ~ Y K [/.-i (Aa) J- (r»P B ) + ^a 2 » +3 / i+1 (Aa) 



a /p« 



3as \r 



?i+i 



+ B W 



(n + !)/»-, (to) A ( r »P B ) _nA»a~»/. +1 (*a) A /2 



and two other similar expressions must vanish, when r = a, for all values of n. 
These three conditions will be satisfied if 



Kfi(ha) = ~-\jj h (ha) 



• • » t « « • 



• • (8): 
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and if for all values of n > 0, 



— ~ ■—A. n f n - l (ha)^B n (n-\-l)f n ^ 1 (ka) = -L n h n ^ n ^(ha) 

- ~A n h 2 a 2 f n+1 (ha) -B n nk 2 a 2 f n+l (ka) = - L n h n h 2 a 2 xfj n+1 (ha) 



>. . . (9). 



J 



These equations (8) and (9) are sufficient to determine the various constants in the 
expression (4) for the scattered sound. In the process of approximating to the values 
of these constants we shall limit ourselves, as before, to the ca.se when ha is a small 
fraction ; in other words, we shall suppose that the radius of the obstacle is small 
compared with the wave-length of the incident sound. We shall also suppose that 
the gas in question is the air of the atmosphere. This will make crvfc 2 a small 
fraction for all wave-lengths within the limits of audibility. 

With this assumption we may write, as before, 

A = £(l-fu*r/c*) (10). 

Eliminating A n between the equations (9), we obtain 

B„{(^+ 1) h 2 a 2 f n+1 (ha)f n ^ (ka) +nf n ^ (ha) k 2 a 2 f n+1 (ka)} 

= -L n h n . h 2 a 2 {xjj n ^ (ha)f n+1 (ha) -y]j n+1 (ha)f n ^ (ha)} . (11). 

Now with the help of a well known result in the theory of Bessel functions it may 
be proved that 

^ (ha)f n+1 (ha)-xj, n+1 (ha)f n ^ (ha) = (2n+ 1) h-**+*>aT<* + *>. 

Hence equation (11) takes the form 

B^(n+ 1)/*V/ W+ ^ = -(2^+ iy^r<* +1 >.a-<* +1 > (12). 

Retaining only the principal term in the coefficient of B w , we may write approxi- 
mately 

whence we have 

**- irh ( n +l)(2w) !"/»_! (ka) ' * h 

to the same degree of approximation. 

Next eliminating B M between the equations (9), we obtain 

— — A M [(w+ 1) h 2 a% +1 (ha)f n - x (ka) + nf n -i (ha) l 3 a%^ (ka)] 

= -Ch n [n^„_! (ha) &V/ M+1 (ka) + (n+ 1) h 2 a 2 xft n+1 (ha)/^ (ka)] . (14). 
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Hence by use of the expressions (13) and (15) § 6 for \fj n (£) andj^(£), we have 
approximately, when ha is small, for all values of n> 0, 

n h 2n+1 a 2n+l 

A n = — t n h n — — p g 3 . — /2 W _i\a * a /^+i (^ a )/fn-i \kci) . . * (15). 

We shall find it necessary to obtain a closer approximation to the value of A x . 
Writing n = 1 in (14) we have 

\A 1 {2h 2 a 2 f 2 (ha)f Q (ka)^f (ha)k 2 a 2 f 2 ^ . (16). 

Now from (15) § 6 we have approximately, since ha is small, 

f Q (ha) = h-'a- 1 (l-iha-ih 2 a 2 ), 2h 2 a% (ha) = 6h~ s a~ s (l+ih 2 a 2 +^hW). 

Substituting these expressions in (1.6) and making use of (13) § 6, we obtain 

2A 1 h~ 3 a~ 3 = ih [l-Zik-^-Zk-^-Wa 2 (^^k~ 2 a~~ 2 ^Lk~ z a~ z ^k"'a~') 

-\h?a? {\i+ 2k~ 1 cT 1 - bJc~ 2 a~ 2 - 6k~ s a~ s + 3iT 4 cT 4 )]. 

Writing Xae~ v * t7r for ka we obtain finally 

A 1 = i^V[V2X- 1 a~ 1 + 3X" 2 a~ 2 +p 2 a 2 (3X~ 2 a~ 2 + 3y2X~ 3 a~ 3 ) 

-i^a 8 (£+ v / 2X- 1 a" 1 + 3^/ 2X- 3 a~ 3 - 3X~ 4 a~ 4 )] 
+ 6|AV[l+|y2X- 1 a~ 1 -|A 2 a 2 (f + 3 x /2X-- 3 a- 3 +3X- 4 a™ 4 ) 

-i^a^^X^a^ + SX^a^-S^X" 8 ^ 8 )] . . . (17). 



By a similar process we obtain from (8), 

Ao = -^ 8 a 8 (l-fA 2 a a +iiA 8 a 8 ) (18). 

§ 8. Having obtained approximate values for the various constants involved in the 
expressions for the secondary waves, we may now proceed to estimate the additional 
loss of energy consequent upon the presence of the obstacle. The method adopted is 
exactly similar to that of which we made use in the case of cylindrical obstacles. 
As above, it is easily seen that the total additional dissipation of energy due to the 
presence of the obstacle is given by 

Ipoffi ^S + i?i9o dS ......... (1) 

where the suffixed letters have the same meaning as in § 4, and the integration is to 
be taken over the surface of a sphere of radius R concentric with the obstacle. As 
before, we shall suppose that R is great compared with the wave-length of the 
incident sound, and yet such that crVR/c 3 is a small fraction. By this assumption we 
are enabled to neglect the imaginary part of KR and also to rega-rd the motion as 
sensibly irrotational at the boundary r = R, 
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Now, at the external boundary r = E, we may write approximately 

Bf = I [A B /r<" +1 > . t "P n (ji) e'^- M >] 

where square brackets are used to denote that the real part only of the expression so 
enclosed is to be taken into the account. 
Hence we obtain 



CO . . co 



R^ = - po(r 2 [A K r <» + V +1 P n (ji) . e ^- m) l Kfc = h t [A^- (a+ V +1 . P M (ji) &«-**>]. 



n-Q n=0 



Combining these results we find 

R(Mi+l>igo)H-/^ • • (2)- 



ft=0 



Again we have from (3) §7, since KR is large, 

K<j&o = X {(2n+ 1) A _1 P n (/x) sin (liR—^nrr) cos (crt + ^nir)}. 



Hence we have 



Rp = ^ oC r I {(2^ + 1) h X P W (/a) sin (/?H— |wr) sin (cr£ + |-n77-)} 



rc=0 



R^o = _ /*, £ {(2n + 1)4 *P n (/x) cos (AK— Jn7r) cos (crt + ^mr)} 



71 = 



Combining the two last results we find 

R(-po<rqo + hp Q ) = po<r S {(-)* (2n + 1) P w (/x) cos (at-hU)}. 



»=0 



Substituting this result in (2) and integrating over the surface of the sphere r = R, 
we obtain 

"(P&o+Poqi) dS = 4ir/> <r i [(-) n A,t w+i r (n+1) cos (<rf-£R) e t( ^ R) ], 

of which the mean value is 

27r Po <r Z [(-Y A n i n+1 h-^~] (3). 

n-0 

This last expression then represents the loss of energy to the primary waves in 
consequence of the presence of the obstacle. From the value of A n obtained in (15), § 7, 
we see that the summation (3) consists of a series of terms arranged in descending 
order of magnitude. Consequently, in determining its value we may limit our 
attention to the first two terms. Hence the total loss of energy to the primary 
waves is given very approximately by 

2Trp crh[A ih~ 2 +A 1 h~ z '] . . . , ■ (4). 
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Now, from (17) and (18), § 7, we have 

[A ^- 2 + Vr 3 ] = iha s {3 v/2 . X-'a- 1 + 6\~ 2 a~ 2 + 3h 2 a 2 (\~ 2 a~ 2 + ^/2\~ z a~ z ) 

+ ^ 3 (|+y2X~ 1 a- 1 + 3 v /2X- 2 a- 2 ~3X- 4 a- 4 )} . . . (5). 

Let us first consider the case when \a is small. Since h 2 /k 2 or ovjc 2 is always small, 
it follows from (5) that in this case we may write approximately 

[A ih~ 2 + M" 3 ] = \a 2 U ^ + 3 v/2o-VYc\ 

When Xa is great, it is necessary to include one other term of (5), and we may in 
general write in this case 

[A ^- 2 +A 1 / i - 3 ] = la 2 (|o- 4 «Vc 4 +3 x /2o-V / yc+6^-) .... (6). 

Comparing this last result with that obtained in the case when Xa is small, we see 
that we may take it as a sufficient approximation in almost all cases. For small 
values of the radius the first term in (6) will be negligible. 

Substituting from (6) in (4) we obtain for the total loss of energy to the primary 
waves the expression 

ip (T 2 /c7ra 2 . (f o-W/c 4 + 3 v /2o- 1 V /2 /c + 6 — 

Now the energy incident upon the obstacle in the primary waves is given by 
|-/) (r 2 /c7ra 2 , and hence the ratio of the lost energy to that incident upon the obstacle is 

3^2(7^/0 + 6— + |crV/c 4 (7). 

The first two terms of this last expression represent the proportion of the incident 
energy lost by friction. The last term of (7) gives the proportion lost by scattering 
to a distance, and is the same as is obtained in the theory of a frictionless air. 

When Xa is small, the most important term of (7) is the second 6 — . Hence we 

ca 

see that in the case of small obstacles the ratio of the lost energy to the incident 
energy varies inversely as the radius of the obstacle, and consequently tends to become 
very great as this radius is diminished. On the other hand, the actual amount of 
energy lost varies directly as the radius of the obstacle, and diminishes with the 
radius. It is to be noticed that in the case of sufficiently small obstacles the energy 
lost to the primary waves is independent of the wave-length of the incident sound. 

When Xa is great, the most important term of (7) is the first 3 A /2cr 1/ V /2 /c. Hence 
we see that in this case the ratio of the lost energy to that incident upon the obstacle 
is very nearly independent of the radius of the obstacle, provided the order of 
magnitude of this ratio is altered by the viscosity. Consequently for sufficiently large 
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obstacles the loss of energy to the primary waves is proportional to the surface of the 
obstructing sphere. There is, it will be noticed, a distinct similarity between these 
results and those obtained above in the case of cylindrical obstacles. 

The expression (7) has been evaluated in a number of different cases, and the results 
are arranged on p. 266 in tabular form. K denotes the expression (7) or the ratio of 
the lost energy to that incident upon the obstacle, and X represents the wave-length 
(measured in centimetres) of the incident sound. 

§ 9. Application of the above to the Problem of a Large Number of Spherical 
Obstacles, — -Let us consider now the loss of energy to the primary waves when these 
are incident upon a large number of spherical obstacles. We shall suppose that there 
are n small spheres per c.cm. ; the validity of our argument will depend on the 
volume occupied by the obstacles being small compared with the total volume. 
Consequently ^nrra s must be a small fraction. 

At a distance r from the centre of any one of these spherical particles, great 

compared with the wave-length of the incident sound, the secondary waves due to 

that particle will be sensibly irrotational, a ad will be given very approximately in all 

cases by 

fa = A / (hr) + A x f (hr) r/x. 

Since hr is great, we may write 

p — ihr p—ihr 

and fa takes the form 

fa = (Aoh^ + AJi-'L^e-^/r; 

which, along the course of the primary waves (/x = —1), reduces to 

<t> 1 = (A&- 1 -A 1 hr\)-e-' kr lr (1). 

Consider now the spheres which occupy a thin stratum dx perpendicular to the 
course of the primary waves. Let P be any point in this stratum, and let O be the 
point where the vibration is to be estimated at a great distance from the stratum. # 

If AP = 2, the element of volume is 2irx.dx.dz, and consequently the number 
of spherical particles in it is 2irnz.dx.dz. Also, if OP = r, AO = — x y then 
r 2 = x 2 +z 2 and rdr = zdz. 

Now by (1) the resultant at of all the secondary vibrations which issue from the 
stratum is given by 



/*00 



2nn dx (Ao/T 1 - AJfh) e~ lhr dr. 

J —x 

Remembering that the angle AOP is to be regarded as very small, we see that the 

* See figure, p. 254. 
VOL. OCX. — A. 2 M 
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resultant at of all the secondary vibrations coming from the stratum dx is given by 

-2im dx (iA*h~ 2 + L x h-*) e«^ . ....... (2), 

where the time factor e l<Tt has been restored. 

Now, from the results obtained in § 7 for A x and A , we find 

tAoA-'+Aifc- 8 = i^ 3 {3 v /2\- 1 a- 1 + 6X" 2 a- 2 +3AV(X- 2 a~ 2 +y2X~ 3 a" 3 ) 



1 

2 



A 8 a 8 ( v / 2X- 1 a~ 1 + 5\- 2 a- 2 ~3 v / 2\- 3 a- 8 )}. 



Except for very minute obstacles, it will be sufficient to write 

lAo^+AA" 3 = Pa 3 {3 v /2X- 1 a- 1 + 6r 2 a- 2 +^V} + £i>>a 8 (£ + f v /2r 1 a" 1 ). 

Substituting this last expression in (2), we obtain for the resultant of all the 
secondary vibrations coming from the stratum dx 

-\mr dxo-a d /c [_(3 JS /2\' 1 a' l + 6\~ 2 a~ 2 +^ 3 a 3 ) + 1 (f + S^X" 1 ^" 1 )] e l( ** +a *>, 
of which the real part is 

—^nwdx.aa z /c{(3^2\~ 1 a^ 1 -{-6X~ 2 a~' 2 -{-^h 3 a d ) cos (hx + crt) 

-(f +3 x /2\" 1 a~ 1 ) sin (Aa? + <rf)} . (3). 

To this is to be added the corresponding expression for the primary wave 

(j> = cos (hx+crt). 

The coefficient of cos (hx + crt) is thus altered by the obstacles in the layer dx from 
unity to 



1 -hnra 2 \ 6 — + 3 J 2o- 1/ V /2 /c + icrW/c 4 
[ ca i j i 



> dx. 



Thus, if E be the energy in the incident wave, we have 

<£E/E = -nrra'U— +3 v /2o- 1/ V / 7c + |o" 4 a7c 4 l dx. 

Integrating this, we obtain 

E = E e~ ax , 

where E is the energy in the primary waves at incidence, and a is given by 



a = nwa 2 U~ + 3^2ct 1/ V /2 /c+z<tW/c 4 ) (4). 



2 M 2 
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If the radius of each obstacle is measured in centimetres, then the reciprocal of a, 
as determined by (4), will give the distance travelled by the sound before its intensity 
is diminished in the ratio of l/e. If the radius of each small sphere is 10~ 3 cm., and 
there are 10 6 per c.cm., then ^mra? will be a small fraction, and the formula (4) will be 
applicable. With these numerical values we obtain, in the case of sound of wave- 
length 50 cm., a = 8*5 10~ 2 . Consequently aT 1 — 11*8 cm. ; hence, after passing 
through a thickness of less than 12 cm. of such a medium, the intensity of the 
sound will be diminished in the ratio of l/e. 

The formula (4) should be applicable to fogs, as we may regard the water particles 
as approximately fixed, since their inertia is so much greater than that of the 
surrounding air. I am indebted to Prof. Lamb for the following information from 
Hann's ' Meteorologie ' : " In a dense fog the amount of water may vary from about 
3 to 10 gr. per cubic metre. Assuming that the diameter of the drops is *02 mm., 
and a cubic metre contains 4*5 gr. of water, this is calculated to give 10 9 drops per 
cubic metre, and therefore 10 3 per cubic centimetre/' With these numerical data the 
formula (4) gives a" 1 =1180 metres, and consequently it follows that the fog would 
not interfere appreciably with the propagation of sound. But if the diameter of the 
drops could be as small as '002 mm., a fog of the same density would contain 10 6 
drops per cubic centimetre, and a" 1 would be nearly equal to 1^ metres, and 
consequently the sound would be damped very quickly by the fog. # On the other 
hand, Tyndall's observations appear to show that the presence of fog is not 
prejudicial to the audibility of sound, f 

The coefficient of sin (hx + o-t) in (3) gives the refractivity of the medium as 
modified by the spherical particles. If S be the retardation due to the spheres of 

the stratum dx, 

8 = ^mr.dx.a? (1 + 3 ^/2 X" 1 ^" 1 ). 

Hence, if /x be the refractive index of the medium as modified by the particles, 

p- 1 = f p (|+ 3 ^/2 X" 1 ^- 1 ), 

where p denotes the ratio, assumed small, of the volume occupied by the particles to 

the total volume. 

Hence finally we have 

lx - 1 =p{i + ^ x /2( v /af/a} (5). 

For sound of wave-length 50 cm. incident upon a medium in which there are 10 6 
spherical particles per cubic centimetre, each of radius 10" 3 cm., we obtain 

^-1 = 3.7. 10" 2 . 

* See, however, note at end. 

f EAYkEiGH, * Treatise on Sound,' Vol. II., p, 137. 
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Note, April 12th, 1910. 

Prof. Labmor has kindly pointed out to me that it is not legitimate to apply the 
formula (4) to fogs without further consideration. Although the inertia of the water 
particles is so much greater than that of the surrounding air, yet in consequence of 
the viscosity of the air it does not follow that we may regard the water particles as 
approximately fixed. I have investigated the problem of a free spherical obstacle. 
The analysis is very similar to that in the problem of the fixed obstacle. The 
secondary waves diverging from the obstacle are affected only in the terms containing 
spherical harmonics of the first order. If \Je t<Tt is the velocity of the obstacle along 
the axis of x we obtain 



■g- Aj/o {ha) + B x 2f (ha) = — ih\fj (ha) + U . 



(1): 



together with 



^A 1 h 2 a 2 f 2 (ha)-B 1 k 2 a 2 f 2 (ka) = -ihh 2 a 2 xjj 2 (ha) ..... (2). 



MtcrU 



j j 



fS#*'jt* \A/h^J m 



■ (3) 



where p rx is the component of the stress across the surface of the obstacle in the 
direction of the axis of x and the integration is taken over the surface of the obstacle. 
The last equation reduces to 

U = £2[U + 3B l {3/ 1 (*a)-/o(*a)}] . ...... (4) 

Pi 

where p x is the density of the obstacle. 

Eliminating A 1 and B x between the equations (1), (2), and (4) we obtain approxi- 
mately when ha is small 



Ull-|&i-V 8 
Pi 



Pi 



Hence, if L be the ratio of the amplitude of the motion of the obstacle to that of 
the waves of sound, we have 

L= {l+AWpo)V 4 a- 4 }^ (5). 

Again, eliminating B x and U from (l), (2), and (4) we obtain approximately 



Aj, = —^Lh A a d k 2 a 2 f 2 (ha) 



pi 



Hence we obtain without difficulty when ha is small 

[A 1 ] = [A 3 ] {l4-^(p /pi) 2 X-%"" 4 }- 1 . ...... (6), 

where [AJ, [A x ] denote the real part of the value of A in the case of the fixed and 
free obstacle. 
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The last relation may be written in the form 



[A,] = [A 1 ] (1-L») . 



(*)• 



It follows that the ratio of the lost energy to that incident upon the obstacle is 
given by 

...... (8). 



i cr V/c 4 + ( 3 v / 2(7 1/2 ^ 1/ 7c + 6 — ) ( 1 -~ L 2 ) 



Extending this result to the case of a number of free spherical obstacles we obtain, 

instead of (4), § 9 

a = ri7ra 2 K. 

where K represents the ratio of the lost to the incident energy given in (8). 

It follows from this investigation that the results obtained in this paper are only 
applicable to fogs if •§■ (po/pi) A.~ 2 a~ 2 is a small fraction. This condition is satisfied for 
obstacles of radius 10~ 2 cm., and also for obstacles of radius 10~ 3 cm. when the wave- 
length of the incident sound is not too long. In the case of obstacles of radius 
10~ 4 cm., however, this condition is no longer satisfied; L approaches close to unity 
for all wave-lengths, and consequently a and K are very small. Hence, if the 
diameter of the drops of water in a fog is as small as '002 mm., such a fog does not 
interfere appreciably with the propagation of sound, and a result is obtained in 
agreement with Tynd all's observations. 

I append a table giving the values of L in a number of different cases. When the 
wave-length of the sound is very great, or when the obstacle is extremely minute, the 
obstacle vibrates with the air surrounding it. 
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